Toy models of crossed Andreev reflection 



m 
O 
o 



R. Melin*W, H. JiraxiW and S. Peysson^ 
^Centre de Recherches sur les Tres Basses Temperatures (CRTBTf 
CNRS BP 166X, 38042 Grenoble Cedex, France 
Instituut voor Theoretische Fysica, Universiteit van Amsterdam, 
Valckenierstraat 65, 1018XE Amsterdam, The Netherlands 
(February 1, 2008) 

We propose toy models of crossed Andreev reflection in multiterminal hybrid structures containing 
out-of-equilibrium conductors. We apply the description to two possible experiments: (i) to a device 
containing a large quantum dot inserted in a crossed Andreev reflection circuit, (ii) To a device 
containing an Aharonov-Bohm loop inserted in a crossed Andreev reflection circuit. 
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I. INTRODUCTION 

Transport of correlated pairs of electrons in solid state 
devices has focussed an important interest recently. Be- 
cause of the progress in the fabrication of nanoscopic 
devices it will be possible to realize in a near future 
transport experiments in which two normal metal or spin 
polarized electrodes are connected to a superconductor 
within a distance smaller than the BCS coherence length. 
Transport theory of multiterminal hybrid structures has 
been investigated recently with various methods (a scat- 
tering approach in Ref. [1], lowest order perturbation in 
Ref. [2], Keldysh formalism in Ref. [3]). Multiterminal 
hybrid structures can be used to manipulate entangled 
pairs of electrons in solid state devices [1-8] and propose 
new tests of quantum mechanics with electrons such as 
EPR [4,5,8] or quantum teleportation experiments [7]. It 
has been suggested recently [9] that proximity effect ex- 
periments at a ferromagnet / superconductor interface 
could be explained by spatially separated Cooper pairs 
in which the spin-up (spin-down) electron propagates in 
a spin-up (spin-down) domain. This gives a strong moti- 
vation for investigating new situations involving spatially 
separated pairs of electrons, which we do in this article. 
One of the possible experiments that we propose involves 
a large quantum dot inserted in a crossed Andreev re- 
flection circuit. This device can be used to probe spin 
accumulation related to crossed Andreev reflection and 
elastic cotunneling. Another device that we propose in- 
volves an Aharonov-Bohm loop inserted in a crossed An- 
dreev reflection circuit. This device can be used to probe 
Aharonov-Bohm oscillations associated to spatially sepa- 
rated pairs of electrons. By crossed Andreev reflection we 
mean the possibility that a spin-up electron from a spin- 
up ferromagnetic electrode is Andreev reflected as a spin- 
down hole in another spin-down ferromagnetic electrode 
on the condition that the distance between the ferromag- 



netic electrodes is sufficiently small. The other possible 
process is elastic cotunneling in which a spin-up electron 
from a spin-up ferromagnetic electrode is transfered as 
a spin-up electron in another ferromagnetic spin-up elec- 
trode. With partially polarized ferromagnet, crossed An- 
dreev reflection can take place in the parallel alignment 
and elastic cotunneling can take place in the antiparallcl 
alignment. 

The theoretical description is based on toy models re- 
lying on a series of simplifying assumptions. The initial 
electrical circuit is replaced by nodes interconnected by 
tunnel matrix elements [10]. Each node corresponds to 
a large quantum dot so that the energy levels form a 
continuum within each node. We suppose that the prop- 
agators within a given node are uniform in space. A 
ferromagnetic node is thus characterized the spin-up and 
spin-down density of state. A superconducting node is 
characterized by the ordinary and anomalous propaga- 
tors. We suppose that the applied voltages are small 
compared to the superconducting gap and that the prop- 
agators at a given node are independent on energy. Each 
node is supposed to be in local equilibrium so that the 
distribution function within each node is represented by 
the Fermi-Dirac distribution. To impose Kirchoff laws we 
determine the spin-up and spin-down chemical potentials 
so that the spin-up and spin-down currents are conserved 
at each node [10]. 
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FIG. 1. The toy model used to discuss the interplay be- 
tween elastic cotunneling and sequential tunneling. The local 
chemical potential at node x is determined in such a way that 
current is conserved. 
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II. M/M/M JUNCTION 

Let us first consider the metal / metal / metal junction 
on Fig. 1 in which we suppose that a large quantum dot 
(node x) is inserted in between two metallic electrodes 
represented by nodes a and j3. We first calculate the 
Green's functions Gij of the connected system that are 
the solution of the chain of Dyson equations given by 
G = <7 + <7(g)£(g>Gina compact notation [11], where S 
is the self energy that contains all couplings of 

W = t a<x C+C x + t x ^ a C+C a + tp !X C^C x + t x ,pC+C0. (1) 

The coupling Hamiltonian (1) is formally equivalent to 
a tunnel Hamiltonian but similarly to Ref. [12] we use 
a non perturbative method based on Ref. [11] that is 
valid also for large interface transparencies. This method 
avoids inconsistencies present in the tunneling Hamil- 
tonian model in the regime of high transparency inter- 
faces [13]. The symbol ® in the Dyson equation in- 
cludes a convolution over time arguments and a sum- 
mation over the labels of the network. Since there is no 
magnetic flux one has t a ^ x = t XM and t Xt p = tp tX . To 
obtain the current flowing from node a to node x we 
evaluate the Keldysh Green's functions given by G + '~ = 
(1 + G R ® E) ® g+<- <g> (1 + S ® G A ) [11,12]. The cur- 
rent is related to the Keldysh Green's function by the 
relation [11] 

I a , x = jj [t a , x G+;-(u) - t x , a G+ : -(u)] du. (2) 
The current flowing from node a to node x is given by 

(3) 

— (M/3-Ma), (4) 
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with V = 1 +ttH x .aPaPx + n 2 t Xt pPpPx- The current flow- 
ing from node x to node (3 is obtained by exchanging the 
labels a and (3 and changing sign. The chemical poten- 
tial p x in the intermediate region is deduced from the 
KirchofT law I a . x = I x ^\ 



Px = 



t a ,xPaPa + t(3, x Pl3Pl3 
t 2 a ,xP<* + t% x P0 



(5) 



As expected if t a>x — tp tX and if p a = p@ we obtain 
Px = \{Pa + Pf))- The conductance G = dI/dV(V = 0) 
is obtained by replacing the value of the chemical poten- 
tial (5) in the expression of the current Eq. (3): 



G = 



e 2 4tt A ti !X tj x p a p pi 



V 2 



1 + 



PaPx + ^ 2 t1 >x Pf3Px 



(6) 

(7) 



We recognize a contribution due to elastic cotunneling 
and a contribution due to sequential tunneling. 




FIG. 2. The toy model used to discuss Aharonov-Bohm 
oscillations. A voltage V\ — is applied on node 1 and a 
voltage V-i — V is applied on node 3. The chemical potentials 
at nodes 2 and 4 are determined in such a way that current 
is conserved. 



III. AHARONOV-BOHM OSCILLATIONS IN 
NORMAL STRUCTURES 



To discuss Aharonov-Bohm oscillations we consider the 
toy model on Fig. 2. We consider a symmetric structure 
in which p = p\ = ps and p' = pi = p4, and we use the 
notation u = ir 2 t 2 pp'. The conductance is given by 



G = 



Au 



h 1 + 4m + 2m 2 (1 -cos$) 



(8) 



which oscillates with the applied magnetic flux. As 
a consequence our toy model can be used to describe 
Aharonov-Bohm oscillations. 



IV. CROSSED ANDREEV REFLECTION 

Now we discuss crossed Andreev reflection [1-3]. We 
consider the same toy model as on Fig. 1 but now node x 
is superconducting. We denote by p\, pi , p^ and pp the 
density of states of spin-up and spin-down electrons in 
electrodes a and (3. We note g the ordinary propagator 
of the superconducting node and / the anomalous prop- 
agator. We show that / and g should satisfy the relation 
f 2 = g 2 . Solving the Dyson equation and evaluating the 
Dyson-Kcldysh equation leads to the transport formula 
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where the determinant T> is given by T> A 
T L ) + n 2 (f 2 - g 2 )T^T L , and where T T 
\t x ,p\ 2 p\ and T l = \t Xta \ 2 p l a + \t x ,p\ 2 pp. The function g l 
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is given by g^ = g + «7r(/ 2 — g 2 )TK The term (9) cor- 
responds to local Andreev reflections, the term (10) cor- 
responds to elastic cotunneling and the term (11) cor- 
responds to crossed Andreev reflection. In microscopic 
models it is possible to show that because of averaging be- 
tween the different conduction channels the conductance 
associated to crossed Andreev reflection is equal to the 
conductance associated to elastic cotunneling [2,3]. This 
cannot be demonstrated in our toy model because we lost 
all information about the spatial dependence of the prop- 
agators in the superconductor. Instead we choose the 
propagators / and g in such a way that the conductance 
associated to elastic cotunneling is identical to the con- 
ductance associated to crossed Andreev reflection. This 
leads to the condition f 2 = g 2 . This condition ensures 
the equality of the Andreev reflection and elastic cotun- 
neling currents also for large interface transparencies. 



V. EFFECT OF AN OUT-OF-EQUILIBRIUM 
CONDUCTOR 

Now we discuss a device containing a large quan- 
tum dot inserted in a crossed Andreev reflection circuit. 
The non equilibrium effects are described by the out- 
of-equilibrium distribution function - l^i,a) where 

[ii^a is equal to the local chemical potential of spin- 
a electrons in node i. The values of p ij(J are deter- 
mined in such a way that current is conserved. The 
geometry of the device is shown on Fig. 3 and the toy 
model is shown on Fig. 4. Without the quantum dot 
the current response is symmetric: for the toy model on 
Fig. 1 with node x being superconducting the crossed 
conductance G a ,p = dI a (V a ,Vp) /dVp is equal to Qp. a — 
dIp(V ,V a )/dV a [2,3]. As it is visible on Fig. 5 such 
symmetry relations are not valid for the device on Fig. 3. 
With the parameters on Fig. 5 we sec that the total den- 
sity of states at nodes 2 and 3 is equal to 1 + y. The 
absolute value of the density of states p can be larger 
than unity. 
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FIG. 3. Geometry of the device in which a large quantum 
dot is inserted in a crossed Andreev reflection circuit. 

Generically if the densities of states increase at one 
point we will be out of the convergence radius of the per- 
turbative series (namely the term of order t 4 becomes 



larger than the term of order t 2 , the term of order t 6 
becomes larger than the term of order t 4 , etc) but the 
values of the tunnel matrix elements are small enough so 
that we are not in this regime. 

v . ©-®-@-@ v j 
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FIG. 4. The toy model corresponding to Fig. 3. Node 1 
corresponds to the spin-down electrode on Fig. 3. Node 3 
corresponds to the spin-up electrode and node 2 corresponds 
to the large quantum dot. 

The limiting cases y = and y = 1 can be understood 
on simple grounds (see Fig. 5). For y = node 1 is a 
half-metal spin-up ferromagnet (meaning that only spin- 
up electrons are present in node 1) and nodes 2 and 3 
are half-metal spin-down ferromagnets. If we suppose 
that V~i = V ^ and V3 = Vs = we see that there 
is no voltage difference between the superconductor and 
node 3. As a consequence the spin-down electrons aris- 
ing from crossed Andreev reflection cannot be transmit- 
ted to node 3 and give rise to charge accumulation at 
node 2. There is no current in this case. If we suppose 
that V\ — 0, V3 = V ^ and Vs = we see that spin- 
down electrons arising from crossed Andreev reflection 
can be transmitted to node 3 because there is a voltage 
difference between the superconductor and node 3. There 
is a finite current, in agreement with Fig. 5. 
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FIG. 5. Comparison between the crossed conductances in 
two situations for the toy model on Fig. 4: (i) Vi = and 
V3 — V (open symbols); and (ii) Vi = V and V3 = (filled 
symbols). The density of states are given by pN = 1, pi,] = 1, 
Pi, I = 0, p 2 ,t = y, P2,i = 1, P3,t = V and p 3>l = 1. The dif- 
ferent curves correspond to different values of the hopping 
parameter: t = 0.11 (□ and ■); t = 0.12 (o and •); t = 0.13 
(A and A): t = 0.14 (V and ▼); t = 0.15 (0 and ♦). 

Another limiting case is y = 1 (see Fig. 5). For y = 1 
node 1 is a half- metal spin- up ferromagnet, nodes 2 and 3 
are normal metal. We find that the crossed conductance 
is finite if Vi = V ^ 0, V 3 = and is zero if Vi = 



3 



and V3 = V 7^ (see Fig. 5). This should be contrasted 
with the case y — where the opposite occurs. The 
case V\ = and V3 = V ^ can be understood from a 
cancellation between the currents due to crossed Andreev 
reflection and elastic cotunneling so that the crossed con- 
ductance is zero unlike the case y = discussed above. 
The case Vi = V^0,V£ = 0is more complex because it 
involves spin accumulation at node 2 in the presence of 
elastic cotunneling and crossed Andreev reflection. Spin- 
up electrons from node 2 are transferred to node 1 be- 
cause of elastic cotunneling. Because of crossed Andreev 
reflection spin-up electrons are transferred to node 1 and 
spin-down electrons are transferred to node 2. The rate of 
the two processes is identical and there is thus no charge 
accumulation like in the case y = discussed above. A 
naive argument would suggest that the spin-up and spin- 
down chemical potentials at node 2 take opposite values 
and that no charge current is flowing from node 2 to 
node 3 but only a spin current is flowing. As seen on 
Fig. 5 this is not the case since we find a finite charge 
current flowing from node 2 to node 3. We suggest that 
the naive reasoning does not work because the distribu- 
tion function at node 2 may be affected by the fact that 
the two electrons arising from crossed Andreev reflection 
have an opposite energy whereas elastic cotunneling is at 
constant energy. This may induce a non trivial distribu- 
tion function at node 2. 
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FIG. 6. Geometry of the Aharonov-Bohm experiment. 
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FIG. 7. The toy model used to discuss the Aharonov-Bohm 
experiment represented on Fig. 6. 



VI. AHARONOV-BOHM EFFECT RELATED TO 
PAIR STATES 

Now we consider the device on Fig. 6 in which an 
Aharonov-Bohm loop is inserted in a crossed Andreev 
reflection circuit and we use the toy model on Fig. 7. 
We first consider that the Aharonov-Bohm loop is a half- 
metal ferromagnet. In this case the only phenomenon 
coming into account is crossed Andreev reflection for 



antiparallel spin orientations and elastic cotunneling for 
parallel spin orientations. We obtain Aharonov-Bohm os- 
cillations with a negative magnetoresistance (see Fig. 8). 
It may be difficult to probe this situations in experiments 
because the phase coherence length in ferromagnetic met- 
als is very small. This is why we consider the case where 
the Aharonov-Bohm loop is made of a normal metal. 
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FIG. 8. Aharonov-Bohm effect related to non separa- 
ble correlations. Node 1 is a half-metal spin-up ferromag- 
net and the Aharonov-Bohm loop is a half-metal spin-down 
ferromagnet: p^f = 1 and pi t ± = 0, i = 2,..., 6. 
We have shown the variation of Gs,i($) — Gg 1 ^ as a 
function of $. The average conductance defined as 
GWa = (Gs,i(* = 0) + G s ,i($ = tt))/2 is positive. The 
tunnel amplitudes are identical for all links of the network: 
t = 0.1 (□); t = 0.09 (■); t = 0.08 (o); t = 0.07 (•). We 
obtain the same oscillations with a parallel spin orientation 
of the ferromagnetic electrodes but with an opposite current. 
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FIG. 9. The same as Fig. 8 but now nodes 2, ..,5 are un- 
polarized: p^-f = p^j = 1 with i — 2, 5. 

Several phenomena come into account: (i) crossed An- 
dreev reflection; (ii) clastic cotunneling; (iii) spin accu- 
mulation in the Aharonov-Bohm loop (the local spin-up 
and spin-down chemical potentials are not equal); (iv) 
electroflux effect similar to Ref. [14] (the local chemical 
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potential oscillates with the magnetic flux applied on the 
Aharonov-Bohm loop). Another possible process occur- 
ring in this structure with highly transparent interfaces 
is that the two electrons of a Cooper pair are transferee! 
into the Aharonov-Bohm loop, couple to the magnetic 
flux and come back in the other electrode. As it is visi- 
ble on Fig. 9 we obtain also Aharonov-Bohm oscillations 
but with a positive magnetoresistance. 



VII. CONCLUSION 

To summarize we have provided a toy model of mul- 
titerminal hybrid structures containing large quantum 
dots. We applied the model to a device in which a large 
quantum dot is inserted in a crossed Andreev reflection 
circuit. We also proposed another possible experiment 
intended to probe an Aharonov-Bohm effect related to 
spatially separated pairs of electrons. These two situ- 
ations may be the object of future experiments. The 
description was based on a toy model in which the initial 
electrical circuit is replaced by a set of nodes intercon- 
nected by tunnel matrix elements. It is expected that 
the qualitative behavior is captured by our toy model. 
The formulation of a microscopic theory is left as an im- 
portant open question. Within a microscopic theory it 
would be possible to discuss rigorously the role played 
by the elastic mean free path and the phase coherence 
length. In our toy model each node is described by a 
single propagator. As a consequence we have lost all 
information about the spatial variation of the propaga- 
tors and we suppose implicitly that the distance between 
the contacts is smaller than the elastic mean free path. 
But qualitatively for the device on Fig. 6 the current is 
proportional to the square of the anomalous propagator 
in the superconductor, and to the square of the ordi- 
nary propagator in the normal metal electrodes. We have 
supposed in our toy model that the distance D between 
the contacts on the superconductor is much smaller than 
the BCS coherence length £bcs and the length R of the 
normal metal electrodes is much smaller than the phase 
coherence length l^. Qualitatively the conductances on 
Figs. 5, 8 and Fig. 9 should thus be multiplied by the ex- 
ponential factors exp (— -D/£bcs) and exp (— R/l^). An- 
other important ingredient is the role played by interface 
transparencies. We could solve the toy model with arbi- 
trary interface transparencies and obtained in some cases 
a crossed conductance larger than e 2 /h with large inter- 



face transparencies. As discussed previously the value 
of the conductance would be much smaller in a micro- 
scopic model because of the exponential dependence of 
the propagators in the superconductor. To gain in real- 
ism it would be also interesting to discuss diffusive models 
in connection with the occurrence of h/e oscillations in 
the conductance predicted from the toy model, h/e os- 
cillations have been also obtained for the proximity effect 
in other geometries (see Ref. [15]). 

The authors acknowledge fruitful discussions with D. 
Fcinberg. 
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